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(57) Abstract 

Methods and apparatus for delivery of ultrasound to the brain, without requiring removal of portions of the skull, call for transmission 
of ultrasound with a plurality of transducers (12) aimed to induce cavitation at least at a selected region of the brain. An excitation source 
(22) is arranged for driving at least selected transducers at differing phases with respect to one another, e.g., to compensate for phase shifts 
(or phase distortions) effected by the skull on the ultrasound output by each transducer. As a result, the ultrasound waves reaching the 
selected region from the transducers arrive substantially in phase with one another, e.g.. within 90 degrees and preferably 45 degrees, and 
still more preferably 20 degrees of one another. 
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Methods and Apparatus for Delivery of Noninvasive Ultrasound 

Brain Therapy Through Intact Skull 



Sponsorship 



The research resulted, at least, in part, from work performed under NCI Research Grant 
No 46627. 
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Reference to Related Applications 

This application claims the benefit of the filing date of and is a continuation-in-part of 
copending, commonly assigned United States Patent Application Serial No. 08/71 1,289, filed 
10 August 21, 1996 (Attorney Docket No. 0092664-0008), the teachings of which are incorporated 
herein by reference. This application also claims the benefit of the filing date of copending 
United States Provisional Application Serial Nos. 60/034,084 (filed 12/23/96), and 60/045,453 
(filed 5/1/97). The teachings of those provisional. appUcations are incorporated herein by 
reference. 

15 

Background of Invention 

The invention pertains to medical systems and, more particularly, to methods and 
apparatus for non-invasive application of focused ultrasound to the brain. The invention can be 
20 used, for example, in the diagnosis and treatment of neural ailments. 

According to the prior art, treatment of tissues lying at specific locations within the skull 
are limited to removal or ablation. While these treatments have proven effective for certain 
localized disorders, such as tumors, they involve delicate, time-consuming procedures that may 
25 result in destruction of otherwise healthy tissues. The treatments are generally not appropriate 
for disorders in which diseased tissue is integrated into healthy tissue, except in instances where 
destruction of the latter will not unduly effect neurologic fijnction. 

■ 

4 

The noninvasive nature of ultrasound surgery has special appeal in the brain where it is 
30 often desirable to destroy or treat deep tissue volumes without disturbing the healthy tissues. 

1 
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Focussed ultrasound beams have been used for noninvasive surgery in many other parts of the 
body. Ultrasound penetrates well through soft tissues and, due to the short v^ravelengths (1.5 mm 
at 1 MHz), it can be focused to spots with dimensions of a few millimeters. By heating 
tumorous or cancerous tissue in the abdomen, for example, it is possible to ablate the diseased 
5 portions without significant damage to surrounding healthy tissue. 

Notwithstanding the potential benefits of ultrasound diagnostics and therapy of the brain, 
it has been commonly accepted that ultrasound cannot be applied through the intact skull. Early 
experiments, for example, showed that ultrasound is strongly attenuated by bone, and that brain 
10 tissue damage close to the skull results from the high temperatures caused by the energy loss. 
Accordingly, all of the ultrasound brain treatments performed so far have required the skull bone 
to be removed prior to the sonication. This makes the procedure invasive, and expensive with an 
added risk of complications. 

15 The requirement of surgical removal of the skull has been the main obstacle that has 

prevented ultrasound therapy to be widely tested in the brain despite the possibility of its clear 
benefits compared to other techniques. Accordingly, an object of this invention is to provide 
improved medical methods and apparatus, for diagnosis and therapy of the brain. A more 
particular object of the invention is to provide improved methods and apparatus for application 

20 of ultrasound to the brain. 

A more particular object of the invention is to provide such methods and apparatus as do 
not require removal of portions of the skull, via craniectomy or other such procedures. 

25 Still another object of the invention is to provide such methods and apparatus as can be 

used to precisely target regions within the brain. 

Still yet another object of the invention is to provide such methods and apparatus as can 
be used to effect heating or other physiologic change at such precisely targeted regions, without 
30 effecting substantial change in the surrounding, or other, regions of the brain or skull. 



2 
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Another object of the invention is to provide such methods and apparatus as can be 
utilized over a wide range of ultrasonic frequencies. 

Still another object of the invention is to provide such methods and apparatus as can be 
implemented utilizing conventional materials. 

Yet still another object of the invention is to provide such methods as can be 
implemented without excessive expense. 



If 
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Summary of the Invention 



25 



The foregoing and other objects are met by the invention, which provides in one aspect 
methods and apparatus for delivery of cavitating ultrasound to the brain, without requiring 
removal of portions of the skull. 



Thus, in one aspect, the invention provides an apparatus for delivering ultrasound, 
through intact skull, to the brain comprising a plurality of transducers and an excitation source 
for driving each to induce cavitation at least at a selected region of the brain. The excitation 
10 source is particularly arranged for driving at least selected transducers at differing phases with 
respect to one another, e.g., to compensate for phase shifts (or phase distortions) effected by the 
skull on the ultrasound output by each transducer. As a result, the ultrasound waves reaching the 
selected region from the.transducers arrive substantially in phase with one another, e.g., within 
90° and, preferably, within 45° and, still more preferably, within 20^ of one another. 

15 

The excitation source drives the transducers to deliver ultrasound to the selected region at 
a frequency ranging from 0.01 MHz to 10 MHz and, preferably, from 0. 1 MHz to 2 MHz. 
Sonication duration for the ultrasound ranges, according to further aspects of the invention, from 
100 nanoseconds to 30 minutes. According to still further aspects, the invention provides for 
20 delivery of ultrasound to the selected region with continuous wave operation or burst mode 
operation, where burst mode repetition varies from 0.01 Hz to 1 MHz. 

Further aspects of the invention provide an apparatus as described above, in which only a 
single transducer is used. 



Still fiirther aspects of the invention provide methods for operating transducer arrays as 
described above. 



These and other aspects of the invention are evident in the drawings and in the text that 
30 follows. 



4 
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Brief Description of the Drawings 

A further understanding of the invention may be attained by reference to the drawings, in 

which: 

Figure 1 depicts an embodiment of the invention and an experimental setup for testing it. 

Figure 2 depicts an embodiment of the invention for application of ultrasound to the brain 
of an animal. 

Figure 3 depicts a phased array for application of ultrasound to the brain in accord with one 
practice of the invention. 

Figures 4A-4H illustrate the ultrasound pressure amplitude distribution in water across the 
focus of a transducer according to the invention at various frequencies, with and without sicull 
sections in front of the transducer. 

Figures 5A and 5B illustrate the effect of applying ultrasound in accordance with the 
invention to brain tissue. 

Figures 6A and 6B illustrate phase errors measured at the focus of ultrasound transducer 
arrays with a piece of skull in front of the transducers. 

Figures 7A-7C illustrate the pressure amplitude profiles across the focus of an ultrasound 
transducer phased array in water, through the bone, and through the bone with phase correction. 

Figure 8 . illustrates the pressure amplitude distribution along the central axis of an ultrasound 
transducer array without and with the phase correction. 
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Figures 9A-9C illustrate the ultrasound pressure amplitude distribution measured across the 
focus of an ultrasound phased array in water, through skull without phase correction, and through 
skull with phase correction. 

Figure 1 0 depicts an embodiment of the invention for delivery of cavitating ultrasound to a 
patient's brain through the skull using a multi-element ^transducer array. 

♦ 

Figure 11 depicts a method for delivery of cavitating ultrasound to a patient's brain through 
the skull using a transducer array. 
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> Detailed Description of the Illustrated Embodiment 

Discussed below are methods and apparatus according to the invention for noninvasive 
delivery of ultrasound through intact skull to the brain. These permit ultrasound propagation 
through skull to eflfect cavitation, without causing undesired heating of the brain or, more generally, 
central nervous system (CNS) tissues. These also deliver adequate ultrasound power to ablate 
tissues, or to otherwise induce changes, at focal points (or regions) within the brain. As used herein, 
"tissue" refers to fluids, tissues or other components on or within a patient's body. 

Figure 1 0 depicts an apparatus according to the invention for delivery of ultrasound to the 
brain. The apparatus 10 includes an array of transducers 12 disposed on or near the external surface 
of the head of a human patient. The array 12 can constitute a single transducer, e.g., a spherically 
curved piezoelectric bowl of the type described below, though preferably, array 12 comprises a 
plurality of transducers arranged in a one-, two- or three-dimensional configuration. 

Referring to Figure 3, for example, in one embodiment of the invention, array 12 comprises 
60 individual piezoelectric ceramic transducers mounted in a bowl of circular cross-section. The 
transducer elements, which can be, for example, 1 cm^ piezoelectric ceramic pieces, are mounted 
in silicone rubber or any other material suitable damping agent for minimizing the mechanical 
coupling therebetween. Transducer arrays of this type are known in the art, as described, for 
example in Fan et al, "Control of the Necrosed Tissue Volume During Noninvasive Ultrasound 
Surgery Using a 16-Element Phased Array," Medical Physics , v. 22, pp. 297 et seq (1995), a copy 
of which is filed as an appendix hereto arid the teachings of which (e.g., at Figure 1 and the 
accompanying text) are incorporated herein by reference. The construction of a spherically curved 
phased array comprising multiple square-element transducers is shown in Figure 1 and the 
accompanying text of that publication. 

In the illustrated embodiment, each transducer of array 12 is independently driven by power 
and the control elements 18-22 to generate ultrasound for transmission through the patient's skull 
into the CNS tissues. More- particularly, the transducers in array 12 are individually coupled, via 
coaxial cables 16, to separate channels of a driving system 18. Each channel of that system 18 
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includes an amplifier and a phase shifter, as shown. A common radio frequency (RF) signal is 
driven to each channel by radio frequency generator 22. Together, the radio frequency generator 
22 and driving system 18 drive the individual transducers of array 12 at the same frequency, but at 
different phases, so as to transmit a focused ultrasound beam through the patient's skull to a selected 

« 

region within the brain. Unlike prior art systems, there is no need to remove portions of the skull 
beneath the array 12, e.g., via craniectomy or other such surgical procedure. 

The radio frequency generator 22 can be of any commercially available type. A preferred 
such generator is available from Stanford Research Systems, Model DS345, The generator is 
operated in a conventional way. so as to generate an excitation signal, which is amplified and phase- 
shifted by the individual channels of driving- system 18, in order to induce the corresponding 
transducers of array 12 to radiate ultrasound (e.g., in the range 0.01 MHz to 10 MHz). 

4 

t 

As illustrated, each channel in the driving system 1 8 includes a radio frequency amplifier. 
These can be any RF amplifiers of the type commercially available in the art. 

The phase shifting component of each channel of driving system 18 is constructed and 
operated in the conventional manner known in the art. Particulariy, each phase shifter shifts the 
phase of an incoming RF excitation signal, received from RF generator 22, by an amount a^, ccj, a^, 
etc., as shown in the drawing. These phase shift factors ttj, a^, a^^ etc., can be pre-stored in the 
channels of driving system 18 or, preferably, generated by a controller 20. That controller 20 can 
be a general purpose, or special purpose, digital data processor programmed in a conventional 
manner in order to generate and apply phase shift factors in accord with the teachings hereof. 

The phase shift factors, a^, ag, etc. serve two purposes. The first is to steer the composite 
ultrasound beam generated by transducer array 12 so that it is focused on a desired region within the 
patient's brain. The component of each phase shift factor associated with steering is computed in 
the manner known in the art for steering phased arrays. See, for example, Buchanan et al, 
"Intracavitary Ultrasound Phased Array System," IEEE Transactions Biomedical Engineering, v. 
41, pp. 1 178-1 187, a copy of which is filed as an appendix hereto and the teachings of which are 
incorporated herein by reference. Array steering, or focusing, is particularly discussed in that article, 

8 
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for example, at pages 1179-1181 and, more particularly, in the section entitled "Focusing 
Techniques," the teachings of which are incorporated herein by reference. 

The second component of each phase shift factor a^, aj, a^, etc., compensates for phase 
distortion effected by the skull in the ultrasound ouput by each transducer. In other words, the 
second component of the phase shift factors compensates for perturbations and distortions 
introduced by the skull, the skin/skull interface, the dura matter/skull interface, and by variations 
in the skull thickness. As those skilled in the art will appreciate, the two components that make up 
the phase shift: factor for each channel of the driving system 18 are summed in order to determine 
the composite phase shift factor for the respective channel. 

The phase corrections that constitute the aforementioned second component of each 
phase shift factor can be determined a number of ways. In one embodiment of the invention, that 
component is determined from measurements of the thickness of the patient's skull under each 
transducer in array 12. Such skull thickness measurements can be made using conventional 
imaging techniques, such as computed tomography (CT) or magnetic resonance imaging (MRI). 

In an alternative embodiment, the aforementioned second component of each phase shift 
factor is determined by placing the array 12 on the patient's head and exciting individual 
transducers with a short ultrasound pulse. The echo back from the inner surfaces of the skull are 
monitored by the transducer array 12. The effect of the skull on ultrasound generated by each 
transducer is determined from those echos in accord with conventionally known relations. 

In still further alternative embodiments, the aforementioned second component of each 
phase shift factor is determined by implanting small hydrophones in the patient's brain. These 
are used to monitor the phase of the ultrasound generated by each transducer, e.g., in a manner 
similar to that described below in connection with Figure 1. 

In lieu of illustrated components 18-22, the transducer array 12 can be driven by a 
driving system of the type disclosed in Buchanan et al, supra e.g. at Figure 2 thereof, the 
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teachings of which are incorporated herein by reference. Such a driving sj'steni would, of 
course, require modification in accord with the teachings hereof in order to incorporate phase 
shift factors a,, ctj, ^3, etc., having first and second components as described herein and above. 

Referring to Figure 1 1, the system 10 is operated as described below in order to deliver 
ultrasound through the patient's skull to induce cavitation at a desired region of the brain. 

In step 24, the transducer array 24 is positioned on the patient's head. This is preferably 
accomplished in the conventional manner known in the art for insuring ultrasound transmission 
to the brain. The array is typically positioned over, and as close to, the region in which 
cavitation is to be induced. However, where intervening or adjacent cranial or CNS tissues 
might be adversely affected, the array can be positioned elsewhere and focused accordingly. 

In step 26, the aforementioned second component of the phase shift factor for each 
transducer is determined. This is accomplished in the manner described above, e.g., by 
individual exciting each element of the array and measuring the echo back. The alternative 
mechanisms described above can also be used to determine those components. Those skilled in 
the art will appreciate that in instances where the alternative mechanisms are used, they need not 
be performed after the array is positioned but, can be performed at some other prior time. 

In step 28, the remaining components of each transducers' phase shift factor are 
determined. Particularly, those components associated with steering the array for delivery of 
ultrasound to the desired region are determined. Such determination is made, as indicated above, 
in the conventional manner known in the art for steering phased arrays. 

I 

. In step 30, the array is excited, e.g., by control and driving elements 18-22, to focus 
ultrasound in the patient's head. As noted throughout, because the invention provides correction 
for phased distortion induced by the skull, that ultrasound can be supplied directly through the 
skull without the need for removal of a piece thereof The ultrasound is applied in doses and 
timing sufficient to induce cavitation in the desired region, which may be, e.g., from 1 mm^ - 1 
cm , or larger. Those skilled in the art will appreciate that ultrasound waves in the frequency 

10 • 
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range of 0.01 MHz to 10 MHz and, preferably, from 0. 10 MHz to 2 MHz can be applied with 
sonication duration ranging from 100 nanoseconds to 30 minutes, with continuous wave or burst 
mode operation. The burst mode repetition varies from 0.01 Hz to 1 MHz. 

In embodiments where the transducer array 12 includes only a single transducer, e.g., a 
10 cm diameter piezoelectric ceramic element as described elsewhere herein, step 26 is not 
utilized. In such an embodiment, the "array" is aimed based on its focal point. This is 
determined as a function of the size, radius of curvature and frequency output of the transducer 
in the manner known in the art. In a preferred embodiment, these factors are adjusted so that the 
transducer can be placed directly on the patient's skull, as above. However, where minor 
corrections are necessary, the transducer can be spaced apart from the skull, as necessary, in 
order to insure proper positioning of the focal point. 

Theory 

When an ultrasound beam propagates to a deep target location in the brain part of the 
energy is reflected back at the skin-skull interface due to the high acoustic mismatch between 
these two tissues. The propagating wave in the skull suffers attenuation losses due to scattering 
and absorption. The acoustic mismatch at the bone-dura interface causes part of the remaining 
wave to reflect back to the skull. The total insertion loss through skull depends on the 
frequency, and can be, on average, about 10 dB, and 20 dB at 0.5, and 1.5 MHz, respectively. 
The wave is further attenuated by absorption (about 5 Np/m/MHz) while it travels through the 
brain to the target volume. 

An ultrasound beam delivered to the brain can effect change in CNS tissues and fluids 
(herein, simply "tissue" or "brain tissue," etc.) by two mechanisms: heating and cavitation. The 
ultrasound beam can heat the tissue temperature due to energy absorption from the wave 

resulting in different degrees of thermal damage to the tissue depending on the temperature 

< 

reached. For exposures of a few secorids, temperatures of above about 60° C are adequate to 
coagulate proteins and thus, necrose the tissue. The induced temperature elevation during short 
ultrasound exposures depends mainly on the absorbed power (<q>) although the shape and size 
of the focal spot can have a significant impact due to thermal conduction. The rate of 

11 
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temperature rise (dT/dt) at the very beginning of an ultrasound pulse can be calculated from the 
pressure amplitude of the field (P), as follows: 

dT/dt = <q>/ PC [1] 

where 

<q>= aPV pv [2] 

p is the density (p^ = 1030 kg/m' , p. = 1380-1810 kg/m') (s= skull, b = brain), 
c is the specific heat of the medium (c^ =3.9 kJ/kg/°C ; c, =2.1-2.7 kJ/kg/°C), 
V is the speed of sound, and 

a is the amplitude absorption coefficient of the tissue (a^= 5Np/m/MHz, a, = 50 
Np/m at 0.5 MHz, a, = 300 Np/m at 1.5 MHz (if all attenuated energy is assumed 
to be absorbed)). 

In order to achieve the same temperature within the target volume as in the skull, the 
ultrasound beam has to be focused to overcome the difference in the acoustic properties. The 
square of the pressure amplitude (P^) is directly proportional to the ultrasound beam area 
allowing the required area gain (AG) to be calculated from Equations [1] and [2] by making the 
rates of temperature rise equal in the skull and the brain: 

AG - (Pb'/P.O = (aX)(pyp.)Xvyv.)(cyc.) = 30 at 0.5 MHz [3] 

60 at 1.5 MHz 

In addition the area gain has to compensate for the energy loss due to the skull and 
attenuation in the brain between the skull and the focal point: 

AG'= {a, e-^P*^)-' = 13 at 0.5 MHz [4] 

= . 250 at 1.5 MHz 



12 
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where, 

is the insertion loss of skull, 
P is the amplitude attenuation coefficient, 
f is the frequency and x is the depth in the brain. 

The total area gain is the product of these two area gains and is approximately 400 and 
15000 at 0.5 MHz and 1.5 MHz, respectively when the focus is located at the depth of 6 cm in 
the brain. These calculations are first order estimates and do not take into account phase shifts 
introduced by the variable thickness of the bone or thermal conduction and perfusion effects. 

The second mechanism, cavitation requires negative pressure amplitudes that are large 
enough to form gas bubbles in the tissue. The pressure wave causes the bubbles to expand and 
then collapse. The collapse of the bubbles causes high temperatures and pressures that can cause 
direct mechanical damage to the tissue. Cavitation can offer more therapeutic options than 
thermal exposures of brain. The cavitation threshold in the soft tissues and in bone appears to be 
similar. 

Thus, focusing a cavitation-inducing ultrasound beam overcomes the attenuation losses 
in the bone and brain, but need not overcome differences in absorption coefficients, as is the case 
with the heat-inducing exposures. The beam area of cavitation-inducing ultrasound propagating 
through the skull has to be about 13 and 250 times larger than the focal area at frequencies of 0,5 
and 1.5 MHz, respectively. These area gains are 30 and 60 times smaller than the gains required 
for induction of thermal effects. The cavitation is not influenced by thermal conduction or 
perfusion effects. Therefore, it is clear that cavitation has significant advantages over the 

■ 

thermal effects. This is particularly true in instances where the ultrasound energy must be 
delivered to small focal regions that require high frequencies. 

Cavitation requires high pressure amplitudes but only short exposure durations, therefore 
cavitational effects can be induced without significant temperature elevation. For CNS tissues, 
sonications with durations of only 1 ms are adequate for bubble formation. The required peak 
intensities at 0.936 MHz during these sonications are measured to be around 4000 Wcm** and 

13 
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2000 Wcm'^at 1 ms and 1 s exposures, respectively. Using these intensity values and the 
average ultrasound attenuation in brain of 5 Np/m and in bone 120 Np/m at 1 MHz, the 
maximum peak temperature elevation in the brain can be estimated (from equations 1 and 2) to 
be about 60** C and 0.1** C during the 1 s and 1 ms exposures, respectively. The corresponding 
temperature elevations in bone, if the focus was in the bone, are 1800° C and 3.6° C. During the 
sonications the temperature elevation in the bone would be reduced proportionally with the area 
gain. These values are frequency dependent. For example, bone heating would be about 13° C 
for 1 ms pulse at 1.5 MHz. This short thermal exposure is below the threshold for tissue 
damage. Thermal exposures can be further reduced using multiple pulses that can be repeated at 
a low frequency (for example 0. 1 Hz) thus, eliminating a temperature build up. 

Examples 

In one exemplary embodiment of the invention, which is illustrated in Figure 1, 
ultrasound was generated using a one-element transducer array having a spherically curved 10 
cm diameter piezoelectric ceramic (PZT4) bowl mounted in a plastic holder using silicon rubber. 
The ceramic had silver or gold electrodes both on the front and back surface. The electrodes 
were attached to a coaxial cable that was connected to a LC matching network that matched the 
electrical impedance of the transducer and the cable to the RF amplifier output impedance of 50 
ohm and zero phase. The matching circuit was connected to an RF- amplifier (both ENI A240L 
and A500 were used in the tests). The RF signal was generated by a signal generator (Stanford 
Research Systems, Model DS345). 

The ultrasound pressure wave distributions were measured using needle hydrophones 
(spot diameter 0.5 and 1 mm) and an amplifier (Precision Acoustics Ltd). The amplified signal 
was measured and stored by a oscilloscope (Tektronix, model 243 IL ). The hydrophone was 
moved by stepper motors in three dimensions under computer control. The pressure amplitudes 
measured by the oscilloscope were stored by the computer for each location. 

A piece of human skull (top part of the head: front to back 18 cm and maximum width 12 
cm) was obtained and fixed in formaldehyde. The acoustic properties of formaldehyde fixed 

14 
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skull and a fresh skull are almost identical. The ultrasound transducer under test was positioned 
in a water tank the walls and bottom of which were covered by rubber mats to reduce ultrasound 
reflections. The tank was filled with degassed deionized water. The hydrophone was connected 
to the scanning frame, and positioned at the focus of the ultrasound field. 

Utilizing the setup shown in Figure 1, the embodiment was tested at four different 
ultrasound fi*equencies: 0.246 MHz, 0.559 MHz, 1 MHz, and 1.68 MHz. The maximum peak 
pressure amplitudes achievable through the skull at the focus of the transducer was measured at 
each fi-equency. A shock wave hydrophone (Sonic Technologies Inc, ) was positioned at the 
acoustic focus. Bursts of 10-20 cycles were used to separate the acoustic signal from the 
electrical interference that was picked up by the hydrophone during sonication. Resuhs of the 
testing are shown in Figures 4A-4H. 

Particularly, Figure 4A illustrates the ultrasound pressure amplitude distribution in water 
at the focal point of the single transducer driven at 0.246 MHz, without the skull section in place. 
Figure 4B illustrates this same distribution when the skull was positioned in front of the 
transducer as illustrated in Figure 1. Figures 4C and 4D illustrate the same distributions (i.e., 
with and without the skull section in place) for a frequency of 0.559 MHz. Figures 4E and 4F 
illustrate the same distributions for a frequency of 1 MHz. Figures 4G and 4H illustrate the same 
distributions for a fi-equency of 1 .68 MHz. 

To demonstrate that the single-element array of this embodiment delivers sufficient 
energy to induce tissue damage through the skull, in vivo rabbit experiments were performed. In 
these experiments a window of about 15x15 mm was created in the top of the skull. The skin 
was placed over the skull opening and the animal was allowed to recover. A minimum of two 
weeks after the surgery the animal was anesthetized again and placed on top of a sonication tank 
as illustrated in figure 2 and the 0.556 MHz transducer was aimed such than the focus was 
located at 10-15.mm in the brain. The skull piece was positioned in between the traducer and the 
animal. A thermocouple probe (0.05 mm constantan and copper wires were soldered together at 
the tip) was placed on the skull bone (on the side of the transducer that is expected to be the 
hottest location) under a thin layer of connective tissue that was still attached on the skull. Then 
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10 sonications at the maximum power level for the duration of 0.2 s were repeated with the rate 
of 1 Hz. The animal position was moved and the sonicatlon repeated four times in the same 
location with a delay of about 5 min between the sonications to allow the bone temperature to 
return to the baseline. During the 10 s of pulsed sonication the bone temperature increased from 
the baseline of about 30*^ C to maximum of 43° C with rapid decay. After the sonications the 
rabbit was taken to a MRI scanner and Tl, T2 and contrast enhanced scan were performed. 
After the imaging the animal sacrificed. 

Figure 5 A is a scan of the rabbit brain illustrating the effect of 10 sonications for the 
duration of 0.2 seconds, with a pressure amplitude of 8 MPa, repeated at a rate of 1. Hz. During 
the 10 seconds of pulse sonication, the bone temperature in the rabbit skull peaked from a 
baseline of about 30° C to a maximum of 43" C, with rapid decay. The Figure is a T2-.weighted 
fast spin echo image across the brain. The arrow in the Figure shows tissue damage at the focal 
point of the transducer. The skull window on the top of the head is facing down and, thus, the 
ultrasound beam propagated from bottom up. Figure 4B is identical to Figure 4A, except insofar 

* 

as it shows the results where the above sonication was repeated four times. 

Further embodiments of the invention utilize multi-element phase arrays of the types 
illustrated in Figures 3 and 10, in lieu of a single transducer. By controlling the phase of the 
ultrasound wave as a function of transducer location, these embodiments eliminate the phase 
distortion caused by the skull and thus, allow accurate aiming and use of higher frequencies, 
thus, pennitting application of ultrasound to induce cavitation through the intact skull in regions 
of 1 mm to 1 cm^. 

Two phased arrays comprising these further embodiments hkd similar structure and the 
same driving hardware; the resonant frequency being their only significant difference. The two 
arrays operated at 0.6 MHz and 1.58 MHz. The radius of curvature of both of the transducers 
was 10 cm and both of them were cut into approximately 1 cm^ square elements, as shown in 
Figure 3. The total number of elements in both arrays was 64 although only 60 were driven in 
the experiments due to hardware limitations. The ceramic bowl was cut using a diamond wire 
saw so that the elements were completely separated by a 0.3-0.5 mm kerf The kerf was filled 
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with silicone rubber that kept the array elements together and isolated them acoustically. The 
silicone rubber allowed the transducer elements to vibrate with minimum amount of clamping. 
Each transducer element was connected to a coaxial cable and a matching circuit that was 
individually tuned. The arrays were similar to the one described in Fan et al, supra, at Figure 1 
and the accompanying text, the teachings of which are incorporated herein by reference. The 
array was driven by an in-house manufactured 64 channel driving system that included an RF 
amplifier and phase shifter for each channel. The phase and amplitude of the driving signal of 
each channel was under computer control, as described in Buchanan et al, supra, e.g., at Figure 2 
and the accompanying text, the teachings of which are incorporated herein by reference. 

In addition to phased arrays configured as described herein and shown in the 
accompanying drawings, phased arrays can also be constructed in accord with the arrangements 
described and shown in co-pending, commonly assigned patent application 08/747,033, filed 
November 8, 1996, the teachings of which are incorporated herein by reference. 

Using the methodologies and apparatus described above, it is possible to produce a sharp 
focus through the skull with the single element transducer when the operation frequency was 1 
MHz or less. The beams had secondary peaks introduced by the skull but the main peak was the 
largest. The location of the peak was shifted by the skull by 1-3 mm from the geometric focus, as 
shown in Figures 4A - 4H. However, the focus was completely obliterated with an operating . 
fi*equency of 1 .67 MHz, as shown in Figure 4H. 

The maximum pressure amplitudes achieved at the maximum power output of the 
amplifier were fi-equency dependent and are given in the table below. The maximum average 
pressure amplitude at the fi-equency of 0.554 MHz was 8,0+/- 0.6 MPa. 
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Frequency CMHz\ 



Average Pressure 
Amplitude fMPa'> 



P+ 



p. 

fMPa'> 



0.248 



3.8 



4.2 



3.3 



0.559 



8.0 



8.9 



7.1 



1.0 



5.9 



6.9 



4.5 



Figure 5 shows the image across the brain for the first of the sonications and demonstrate 
tissue damage indicated T2 changes. The tissue damage was also visible in Tl images with and 
without contrast enhancement. 

To measure the phase distortion caused by the skull, a hydrophone was placed in the 
geometric focus of the array under test. The skull was placed between the array and hydrophone 
and each transducer element was powered separately in sequence while recording the time 
difference between the reference signal and the acoustic wave at the focus. This was done with 
both of the arrays. The phase changes required to correct all of the waves to arrive at the same 
phase at the focus are plotted in Figure 6. 

To investigate the effect of the phase correction the pressure amplitude distributions were 
measured in water by scanning the needle hydrophone. The main impact of the phase at 0.6 MHz 
was in the location of the focus which could be corrected back to the geometric focus. This is 
shown in Figure 7A, which illustrates the pressure amplitude profile across the focus of the 0.6 
MHz phased array in water. Figure 7B shows the pressure amplitude profile across the focus 
through bone. And, Figure 7C shows the pressure amplitude profile through bone when a phase 
correction according to the invention is used. Figure 8 likewise illustrates the pressure amplitude 
distribution along the central axis of the array with and without phased correction. The 
magnitude was reduced to 33 % and 40 % of its water value without and with the phase 
correction, respectively. The importance of the phase correction was demonstrated more clearly 
with the higher frequency array. With this array the focus was completely destroyed by the skull 
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(Figure 9b). However, when phase correction was introduced, the focal spot was returned into 
its original shape (Figure 9c) with the half-width of the focus of about 1 mm. The insertion of the 
skull reduced the peak pressure at the focus of the 1.58 MHz array to about 5 % of its water 

ft 

value when phase correction was applied. 

The embodiments of the invention discussed above and shown in the drawings provide 
improved methods and apparatus for neural diagnosis and therapy through application of short, 
high intensity ultrasound beams that induce cavitation at selected locations within the brain. 
These and other embodiments can be beneficially used to deliver focused ultrasound beams to 
the CNS tissues and fluids, thereby, permitting their ablation or other physiological modification. 

f 

Thus, for example, the embodiments can be used to ablate tumors, cancers and other undesirable 
tissues in the brain. They can also be used, for example, in connection with the technologies 
disclosed in copending, commonly-assigned U.S. Patent Application No. 08/711,289 (the 
teachings of which are incorporated herein by reference) for modification of the blood-brain 
barrier, e.g., to introduce therapeutic compounds into the brain. Because they do not require that 
portions of the skull be removed, the embodiments permit the foregoing to be performed 
noninvasively. 

The results also show that adequate ultrasound transmission can be induced through 
human skull to induce cavitation in vivo. This can be done with single element applicators, e.g., 
preferably at frequencies less than 1 MHz and at higher frequencies with phased arrays that 
correct the phase distortion caused by the variable thickness of the skull. 

The maximum pressure amplitude of 8 MPa induced through the skull at 0.559 MHz was 
able to induce cavitation damage in vivo rabbit brain. This value was reached through an area of 
10 cm in diameter to a focal spot diameter of about 5 nmi (50 % beam diameter). If the whole 
available skull surface around the brain is utilized, then a window of at least three times larger 
could be used. In addition, the geometric gain would allow the peak power through the skull to 
be increased. Acoustic power up to 30- 80 W/cm' of the transducer surface area for continuous 
wave sonication can be generated by ceramic transducers. Higher peak powers could be achieved 
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with the pulsed sonication used for induction of cavitation. Thus, it is estimated that much higher 
pressure amplitudes than measured here can be induced in the brain through the skull. 

The values measured in connection with the foregoing compare favorably with the 4 MPa 
5 that was reported to be the threshold value in vivo muscle at 0.6 MHz and a value of 8.5 MPa at 
0.936 MHz in vivo rabbit brain (the threshold at 0.6 MHz would be lower since it has been 
shown to decrease with frequency). Thus, the results demonstrate that adequate ultrasound 
transmission through skull can be generated to induce cavitation in the brain. 

4 

10 Our results demonstrate that low frequency beams can be focused through the skull, 

though, the focus may be shifted from its geometric location. Therefore, it can be helpful to 
detect the focal spot location in the brain prior to the therapy exposure, e.g., using magnetic 
resonance imaging to detect the local temperature elevation or cavitation in the brain at exposure 
conditions that are below the tissue damage threshold. For example, low power test exposures 

15 can be delivered through the skull while using MRI to detect the location of the focal spot. 
Based on the imaging information the location can be corrected to overlap the target volume 
prior to the therapeutic exposure. 

The phase measurements with the arrays support the observation made with the single 
20 element transducers shovwng that at 0.6 MHz 80 % of the phase errors caused by skull are less 
than Poland thus, each wave is adding to the pressure wave at the focus. However, at 1.58 MHz 
over half of the waves had phase shifts that caused the waves to arrive out of phase at the focus. 
This observation can be explained by the difference in wavelength that is 2.50 mm at 0.6 MHz 
and 0.95 mm at 1.58 MHz. 

25 

The possibility of inducing selective thermal damage at the focus, without damaging the 
skin or brain surface, may be possible due to the small focal spots achieved with the phase 
correction. However, the thermal exposures have to be short to reduce blood flow and perfusion 
effects that are strong in brain tissue. The sharp temperature gradients at the focus transport 
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more energy away from the focus than in the bone where the beam is wide and the gradients 
shallow. At 1 .58 MHz, full utilization of the skull surface may provide marginally adequate 
geometric gains to overcome the skull heating problem. However, at lower frequencies 
especially around 0.5 MHz the focal brain tissue thermal therapy seems feasible although not as 
5 likely as utilization of cavitation effects. 

The results demonstrate that the effects of the skull to the beam shape can be eliminated 
using a phased array with proper phase corrections. In the example above, the phase correction 
was calculated from hydrophone measurements. The same corrections can be made by 

10 measuring the skull thickness from CT or MRI scans and then calculating the phase correction 
required for each array element. The same may be accomplished by sending a short ultrasound 
pulse from each or selected elements of the of the phased array and then listening for the echo 
back from the inner surfaces of the skull or other structures in the brain. The effect of the skull 
on the wave propagation at each location could then be calculated. This can also be done before 

15 therapy by mapping the skull effect using ultrasound. 

*• 

Although good results were achieved with only 60 transducer elements in the phased 
array, still more and smaller elements may facilitate moving the focal spot inside of the brain. 
Similarly, at higher frequencies, smaller elements may allow better phase correction further 
20 reducing the losses induced by the skull. 

Thus, the geometric gmn of about 20 that is required to compensate for the losses caused 
by the skull can be easily achieved by focusing. This is larger than the gain of 10 required to 
compensate the average losses. This indicates that adequate power for induction of cavitation 
25 can be delivered using phased arrays through the skull even at frequencies that are too high with 
a single element applicator. 

* 

In summary, the invention provides methods and apparatus for noninvasive diagnosis and 
treatment of the brain using cavitational mechanism and pulsed ultrasound. It permits adequate 
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power transmission through the human skull can be induced to cause tissue damage while 
keeping the exposures in the overlying tissues below the cavitation threshold. 

Those skilled in the art will appreciate that the invention can be applied for purposes of 
5 tissue ablation, as well as in other procedures where focussed ultrasound is desired. These 
include opening the board-brain barrier, activation of therapeutic agents, occlusion of blood 
vessels, disruption of arteriosclerotic plaques and thrombi, etc. It will alsQ be appreciated that 

the invention can be applied for treatment of humans, rabbits and other animals. The 

.f 

embodiments discussed above and shown in the drawings are illustrative only. Other 
10 embodiments, incorporating substitutions, modifications and other changes therein, fall within 
the scope of the invention. These include embodiments with transducer arrays of different sizes, 
shapes and numbers of elements, as well as embodiments with different amplification and 
driving systems. In view of the foregoing, what I claim is: 
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Methods and Apparatus for Delivery of Noninvasive Ultrasound 

Brain Therapy Through Intact Skull 



Buchanan et al, "The Design and Evaluation of an Intracavitary Ultrasound Phased Array for 
Hyperthermia," IEEE Trans. Biomed. Eng. . v. 41, pp. 1178 - 11 87 (1994) 
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Design and Experimental Evaluation of 
an Intracavitary Ultrasound Phased 
Array System for Hyperthermia 

Mark T. Buchanan and Kullcrvo Hynynen 



Abstract — For evaluating ihe feasibilicy of crcailng prostate 
cancer, a 64-element linear ultrasound phased array applicator 
for intracavitary hyperthermia was designed and constructed. 
A (S4-channel ultrasound driving system including amplifiers, 
phase shifiars. and RP power meters was also developed to drive 
(he array. The design of the array and driving equipmcni arc 
presented, as are the results of acoustical field measurements 
and in vino perfused phantom studies performed with the array. 
Several techniques for heatinq realistically sized tumor volumes 
A-ere also investtgaicd, including single focus scanning and two 
techniques for producing multiple stationary foci. The results 
show that the operation of the array correlated closely with 
the theoretical model. When producing a single stationary focus, 
the array was able co Increase ti;isue temperature by ll'C Of 
vitrA in perfused phantom. With some minor improvements in 
array design, intracavitary phased arrays could be evaluated in 
a clinical envlronnittiit. 

L INTRODUCTION 

INTRACAVITARY ultra.<;ound arrays offer un attractive 
means of inducing local hypenhermia in deep-seated tu- 
mors located near body cavities. By locating the radiators as 
dose as possible to the treatmeni site, problems frequently 
encountered wiih external techniques, such as blockage of 
the acoustic window by bone or gas. or simply the inability 
to attain adequate encrey penetration, can be avoided. Early 
results using mulliclement. nonfocuscd arrays of half-cylinder 
transducers operating at 1.6 MHz suggest thai such arrays can 
be clinically useful in the treatment of prostate cancer (1). 

The proximity of the prostate lo the rectum wall mukes it 
a good candidate for heating using intracavitary ultrasound 
radiators. Since the prostate is located very near the anus and 
only tnillimeters away from the rectum wall, an applicator can 
be easily located close to the prostate. The prostate is one of 
the most easily accessible tumor sites, and one thai affects a 
large enough population of patients to be potentially clinically 
useful. As such, most of the experiments were conducted with 
the goal of heating the prostate in mind. 

While the nonfocused arrays have shown con!>iderable po- 
tential in heating the prostate, they have two primary lim- 

Manuicript received November 17. N92: revised Aucusi 9. 1994. This 
work w»$ supponed by gram number ROI Ca 489.^9 Vroni the Nuiional 
Cancer Insliluie. 

The autbori tn with ihc Dcpanment of Radfotoev. Division of MRl 
Brigham and Women'* Hospiul. Han.'ard .\tedlcal School, Boston. MA 021 15 
USA. 

IEEE Ug Number 4406164. 



iiations; the depth to which they can effectively heat, and 
their limited ability lo control the power deposition fields 
[2]. I3|. Both deeper penetration and Increased control over 
the power depo.<itinn pattern can he achieved using linear 
phased arrays. Since pha.^d arrays can focus their radiated 
energy, they ihcorciically could heal tissues to therapeutic 
tempernturcit deeper ihnn nonfocused arrays. The ultrasound 
pou'cr deposition piiitern can be electronically tailored as 
necessary to produce therapeutic leinperaiures wiihin a desired 
volume. In ;ncas whiirc heaiins would be undesirable, the 
phased ;jrray could lake advantiiiic ol" destructive imcrfercnce 
to minimize power deposition at those locations. The ease 
by which ihe power deposition pattern can be electronically 
uttered in real linK provides a means for the conipensaiing for 
varying physiologicul parumcieis. patient posilioning. and for 
minimi'/.ing patient discomfort. 

Several types of ultrasound phased arrays have been 
proposed or built for hyperthermic purposes. These include 
Umemuni and Cain's sector-vortex and concentric ring appli- 
cator.^ |4J, (5J. the cylindrical section applicator developed by 
Ebbini ct al. (6), 1 7). as well as the tapered array developed 
by Benkescr <r/ al. f8j. Each of these arrays is composed of 
anywhere from 16 to 64 individual elements and operates at 
frequencies between 500 kHz and 750 kHz. While these arrays 
show significant puieniiai, they arc meant to be used in external 
applications and ihercforc are unsuitable for intracavitary use 
in their reported configurations. 

Previously, a theoretical study on the feasibility of intracav- 
itary pha.scd arrays using half-cylinder elements had been done 
by Diederich and Hynynen [31. Based on many acoustical and 
thermal simulations, it was concluded that a practical array 
would be composed of 30 half-cylinder elements with 2.5 
mm center-to-centcr element spacing operating at 500 kHz. 
The sludy predicted that the grating lobes formed with this 
array could be negated by using surface cooling. One of 
the major con.'Jidcrailons in specifying this design was not 
the performance of the array, bur the apparent high cost of 
the amplifiers necessary to drive the system. The desire to 
minimize the number of amplifiers led to the specification of 
a 30-element array. 

fnitial array designs were based on some recommendations 
from the theoretical .study and expanded upon in an attempt 
to construci a practical intracavitary phased array. It became 
apparent that more than 30 amplifiers would be tiecessary for 
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a useful array. This led to ihc developmcm of the 64-chsnncl 
driving system to overcome the limitations imposed by the 
number of amplifiers. This allowed for the evoluiion of a 6-i- 
channel phased array w«h 1.73 mm cfnicr-tD-centcr spacins 
and a total array lencth of 1 10.5 mm operaimg ai 500 kHz 

U. .\lATERiUS AND METHODS 

A. Array Consit ucitnn 

A 64-elemeni array was constructed using half-cyhnder 
transducers operating in iheir resonant radial mode at 500 kHi. 
The array was made by slicing washer-shaped elements wiib 
a diamond wire saw (Laser Technologies. North Hollywood. 
CA) from 15 mm O.D. bv 30-mm Ions cvlindcrs of PZT- 
4 material (EDO. Sali Lake City. LT) wjih silver electrodes 
plated on both the inner and outer wall sunaces. The transducer 
slices were glued logciher using a silicone adhesive iDow 
Coming. Midland. .MI) with 0. 17-mm thick silicon rubber 
spacers (SPC Technology. Chicago, ID between each element 
(producins a dead space of 0.23 mm between elements). The 
Stack of elements ^at then cut in half along the axis of the 
cylinder and the two half-cylinder sections glued together to 
form the full auay. The array was bonded to a brass shell 
to foim ihe complete applicator, as shown in Fig. I. V/ires 
were soldered to the inside wall of each array element that 
extended the lensth of the shell to the handle where (hev were 
each connected to a 2-m RG-178 coaitial cable. The'electrodes 
on the outer surface of the array were electrically connected by 
single wires embedded in heads of silver epoxy (Chomencs. 
Wobum. MAi running along both edges of the array.. 

B. Driving Hardware 

One of the primary disadvantages of phased arrays is the 
increased complexity of the driving equipment. Due tu a lack 



of commercially avanabie equipment to drive the array, a 
64-channci ampiiner system was developed. This system is 
compo>ed nf phase shifters/duty cycle controllers, amplifiers, 
RF power mcrcfs. ^r.d impedance motchins networks (sec 
block diacra:r. tn ftz and is controlled by an IBM PC 
comoaiibtc comoutcr . u d dtgtial I.'O intcrf'ace 

Tnc amonncrs arc based on a switching MOSFET design 
ano are oesicncd arcjr.d inicmational Reciilier's IR2110 dual 
.M05FET Orivcr Each of the amplifiers is capable of deliv- 
ering up to 16 W of RF power at 500 kHz mio a 50- W load 
from each cnanncl or a total output power of about 850 W. 
The amplifiers convert digital lojiic input signals mio high 
power Sine wave.s v^hilc preserving the phase of the input 
-iignal The jmpiitudc of the output signal is controlled by 
iwo mclhodi All of ir.e ampiiner outputs can be controlled 
simulianeousiy by adjusiine the output voiiaee on a 1000 W 
DC supply, or ihe ouip'Ji of each channel can be individually 
conifolled by varying the duty cycle of the input signal. The 
result of reducing the duty cycle of the inpui signal (the duty 
cycle IS the percent of "on" time of (he input signal per clock 
cycle) IS a correspondme decrease in the amplitude ox the 
output .signal. 

Since (he anipliners require digiiai input signals, the phase 
shifiins and duty-cycie conuoi is implemented using digital 
;ouniers |lO). (111. These circuits provide 22.5' phase shift 
resoluuon from 0-360". and 16 steps of duty cycle (amplitude) 
that aliows {he ooipui lo van* from 0-100% of the output 
power, as allowed by (he voltage from the DC supply. While 
the phase shift resolution may seem somewhat limited, studies 
have shown that it is sufficiem for the uses described here 
(111. {I2i- 

The actual RF output power from each channel is monitored 
by power meters that measure both the ror*ard and reflected 
RF power 1 13). Tne power meters.. which were also designed 
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and constnict^d in house. aJiou for ver.ficai.on of (he ourpui 
power on each of the empliner channels, and allow for 
the monitonnf of ihc ouipuii .o ihai possiblv danserous 
fault conditions do noi occur. B) monuormg the amouni of 

matchine and the 

inay operation on be e.sily r:i..r.ito«d so ihat p'eak .rrav 
p«rtonnance can be mainuined. 

Finally, lo maximize po«er KiMitr betueen the amplifiers 
and airay elemems. a jimple LC :r.ttd.r.ce m=tchms network 
wa$ used on each cnannel b«t«ttr. oo«er meter's and the 
atiay. This not only maxim.:,, ^.y.,: transfer, but ensures 
greater umtonnity m the acou<;;;i: :c,vjt signals for a given 
excitation phase and impl.tud. -.tr-.c *,perimentallv"by a 
needle hydrophone, Uuho.t c,e:-.-.-L. .-na.ch.ng. rhe' phase 
hift and amplitude distomor, ci..>t: v. ,he imoedatwe dir. 
ferences beiueen elements :a v.t ..Ti> would have to be 



characterized using a hydrophone and the results used to 
compensate the excitation signals. 

C Acnujlic Modelint 

The 3-D ultrasound modeling program used was originally 
developed by Diedcrich and Hynynen (3) and is only briefly 
introduced here. The routine models A* cylindncal rad.aiors of 
nnite length, radius and separation as each of che ultrasound 
transducers in the array, and assumes a target media having 
similar acouMical propcnies as tissue (Tabic I). The routine 
models the suriace of each of che cylindrical elements as an 
evenly spaced grid of simple hemispherical sources and uses 
Huygcns principle to sum the contributions of each source at 
each poim in the field. All of the simulations were done wi(h 
the sources spaced on a 1/32 grid. 
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The acoustical pressure field was calculated in (he z-piane 
using 



11- 



(I) 



where 



^ r. 



(2) 



and n, is ihe number of sources in (he z.direction. Hf is ihe 
number of sources in the theu direction. P, is the pressure 
from a single source (Pa). P^a is the pressure amplitude at the 
surface of the source (Pa), r„ is the radius of the source (m). r 
is Ihe radial disurtce from' the center of the source (m). I is the 
u^avclength (m). 6 is (he phase of the excitation signal (rad). 
/ is the operating frequency (Hz), and q is the attenuation 
coefficient (Np/m). 

D. Focusing Techuiifues 

The single focus ca^ie is the simplest Jorm of focusing that 
can be done with a phased array. The single focus is produced 
by senine ihc phases of the driving signals so thai constructive 
interference occurs at ihe desired fncal position. The phaae of 
each of the driving signals for an array with S elements can be 
calculated from the differences in the path lengths if. between 
each array element and the focal position bv 



CO, = 2ir 



(3) 



where /, ii the phase tradiansl of the 7'" clement. / is the 

wavelength (mj. » = 1. 2 ;V and rri = 0. 1.2 The size 

of the focus produced with this technique is usually too small 
to heat an entire tumor volume, and therefore other techniques 
must be used to heat larcer volumes. 

To heat larger volumes, two techniques were investigated: 
single-focus scanning and multiple focusing. The single focu$ 
scanning routme simply stepped a single focus back and forth 
along a predetermined length of the array. The focal depth >^'as 
kept constant while the focus was stepped every 300 ms to the 
next position in the scan. Delivered power was maximal at the 
edges of the scan but was reduced to 64% 01 maximum power 
at the center of the scan to flatten the temperature distribution 
in the perfused phantom expenmems (the power distnbution 
was CAperimenially determined). 

The other tec'nniquc. multiple focusing, simultaneously pro- 
duces more than one focus within the target volume. The driv- 
ing signals necessary to produce multiple foci were calculated 
using two techniques:* split focusing and the pseudo-inverse. 
To create multiple foci with the split focusing technique, the 
array was divided into subairays. each of which produce a 
single focus m the same manner as previously described. The 
pseudo-inverse method, developed by Ebbini and Cain (15). 
uses a series of control points that represent the magnitude 
of the ultrasound field at given points. A brief summary of 
Ebbini and Cain's technique follows. 
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The technique solves the Rayleigh-Sommerfcld integral for 
an arrav with N elements, and a field of M control poini^ 



.V 



T 



(4) 



where p is the density of the propagation medium, c is the 
speed of sound in the medium, k is the wavenumbcr. S' is the 
surface of the source. Vn is the panicle velocity norma) 10 the 
surface of the source, and Tm- is the distance between j 
control poin( (r^) and the surface cf a source (r„). m = l, 
2» ,„, .V/» and r» s= 1« 2» A'. 

To simplify programming of this algorithm. (4) can be 
described bv 



Hu^p 



(5) 



where u is the excitation vector, p is the complex pressure 
at the control points, and H is the remainder of (4). For 
iniracavuary hypenhermia uses, the number of elements in 
the aa;iv (A*) is alwavs ereater than the number of conuol 
points (.^/). This leads to an undcrdctcrmmcd system of 
equations wjih an mfinite number of solunons. The minimum 
norm soiuiion (u) can be determined by using a least squares 
approximation of f3) 



(6) 



where //" is ihe complex conjugaie of H. With this technique, 
all of ihc array ricmcnis contribute to all of the foci, unlike 
(he sphi focusing technique where each element contributes to 
only one focub. 

III. E.KPERIMENTAL TECHNIQUES 

A Aiousncal Efficiency Measurements 

The acoustical output of the cylindrical transducer elements 
was measured using a radiation force technique (161. The 
arrav was placed in a brass cone with 45* sidca co reflect 
ihe radially emitted ultrasound fields down into the acoustical 
absorber [ 17|. Tnc force on the absorber was measured using 
a Metier .^E 160 (Hightstown, NJ) microbalance, while the 
RF power was measured using a Hewlett Packard 438A Rf 
power meter and □ Werlaione C2625 (Brewster. NY) dual 
directional coupler. The efficiency was calculated as the ratio 
of the acoustical power to the RF electrical input power. 

B. Uhrasound Field Measurements 

The ultrasound fields were mapped in a tank of degassed, 
deionized water by mechanically scanning a thermocouple 
embedded in a small (2-mm diameter) plastic sphcit. The ther- 
mocouple was positioned by a three-axis computer CDntrollcd 
scanning table. The applicator was mounted on a rotational 
device that allowed measurements to be made in a radial arc 
around the array by rotating the array. Measurements were 
made on a 1 x 1 mm grid with the recorded data being the 
average of three consecutive measurements. 
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C. /n Vifro Kidney Experiments 

Alcohol fixed canine kidneys were used as phanioras for 
studying the heating characierisiics of the airay. The kidneys 
had previously been prepared as dcjcribtd by Holmes er 
ai 1 18), and were rehydraicd prior to use The cxpcrimcnis 
were conducted ai room icmpcraturc using degassed, dcionizcd 
water as the perfusate. A meiering pump (Fluid-Mctcrinx Inc.. 
RHICKC, Oyster Bay, NY) connected !o the renal aneiy 
circulated water through the kidney while ihe renal vein was 
allowed to drain into the unk. The kidney was held in place by 
gently sandwiching it between two PVC membranes mounted 
to a Plcxielas frame. The applicator was firmly clamped to the 
frame to maintain a fixed distance bciwccn the surfaces of the 
kidney and array. Fig. 3 shows a diagram of the experimental 
setup. 

Temperatures were measured in the kidney using either 
seven sensor probes sutured in place perpendicular to the 
anray through the focal region, or by one or two jmgle sensor 
thennocouples pulled along a path parallel to the array.. In 
experiments conducted with the multiple sensor probes, the 
kidney was exposed to uhrasound for a total of 10 min with 
temperature mcasurcmcntj occurring evcr>' 30 s. During each 
temperature measurement, power to the array was disrupted 
for approximately four seconds; (one second prior to the 
first reading, and up lo three seconds to read all of the 
thermocouples). Since this technique docs not give very nne 
spatial detail, of the temperature distributions, the putt* back 
technique was more frequently used. 

The pulUback experiments were conducted by pulling one 
or rwo single uncoared thermocouples (0.05-mm wirej by a 
computer controlled stepper motor along a track parallel to 
the array through the kidney tissue. The temperatures were 
measured every I mm and the average of three readings was 
recorded. Prior to each of the experiments, a baseline icm-' 
perature was established along the path of the therTnocouplcs. 
The kidney was exposed to ultrasound for 20 min to allow 
the temperatures to reach steady-state before the icmpcrature 
profiles were measured. The difference between the baseline 
and final measurement was used to calculate the temperature 
rise. The kidney was allowed to cool to room temperature 
(typically, 30 min) before to the next experiment w»s .itarted. 

The thermocouples were mainly located in the medulla of 
the kidney. The steady-state temperatures In the medulla have 
been shown to be a strong function of the flow into the kidney 
(19]. The flow values were kept relatively low m order to 
simulate the perfusion in the prostate (20], [21 j. 

IV. RESfLTS 

A. Acoustical Efficiency Measurements 

The acoustical efficiency of the 500-kHz PZT-4 transducer 
material was measured for a full cylinri«r. 30 mm lone, and 
for two half-cylinder arrays with 2.5 mm and i.8-mm center- 
lo-ccntcr element spacing (2.2 mm and 1.5 mm clement 
widths) with total array lengths of 24.7 mm and mm. 
respectively. Small arrays were tested so that the resulting 
field would be collimaicd. unlike a sinele arrav element 
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I width < wavelcngihi. which would produce nearly spherical 
*ravcfronii. The colhmation is necessary to assure the waves 
reflected by the rcfiecior are normally incident to the surface 
of the absorber 

The results of these experiments arc shown m Fig. 4. None 
of the arrays exhibited very high efnciencies despite the low 
operating frequency Though not shown here, some results 
were verified using caiorimetric techniques. While the 30-mm 
long full cylinder element had an efficiency of 7]^. the arrays 
exhibited efficiencies of 27% and 1 1% for the 2.2-mm and 1 .5- 
mm arrays, rcspecitvciy. These efficiencies indicate mat about 
SO^c of ihe eleciricai power delivered to the 1.5-mm wide 
elements used in the pha.scd arrays is lost either as heat m the 
transducer or in the electrical matching and transmission lines 
between the amplifier and the array. 

fl. Uhrasound Field Measurements 

The effects of element spacing were .studied by measuring 
(he ultrasound f^cld dis(rib<itions produced by two lo-element 
arrays wnh l.ifk-mm and 2.5-mm cenier.ro-cemcr clement 
spacing Pig. 5(3) and ibi shows the acoustical field plots made 
in a plane parallel to the array and normal to the surface of 
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the half cylinder eiemcnis for two 16-eleTnent arrayb with I.S- 
mm and 2.5-mm element spacing along with their rcspcciive 
simulated fields (Fig. 5(c) and (d)). Both array? were focused 
30 mm from the array sunacc and 15 mm above the ceniral 
axis of the array. Tne most significant effect illustrated here is 
ihe generation of the grating lobe by the l6-e)emcnt array with 
2.5-mm ccnter-io-cenicr clement spacing (Fig. 5tb) and (d)). 
Since such a large grating lobe not only reduces the power 
of (he focus but also could couse heating in unwanted areas, 
arrays must be designed lo mmimize grating lobe formation. 
The array with KS-mm cenier-io-cenier spacing produced 
much better field disrhbutions than the array with 2.5-mm 
cemer-to-center spacing. Therefore, in order to have center- 
lo-cenier spacing I.S mm or smaller. 1.5-mm wide elements 
were used in (he final array design and the dead space between 
elements was reduced to 0.23 mm. 

The complete measured and simulated plots of the ultra- 
sound field produced by the 64-clemcm array with 1.73-mm 
element spacing are shown in Fig. 6(a) and (b), respectively. 
The array was focused 40 nrun deep, and 30 mm from the 
central axis. Note that only a small grating lobe is produced 
by this anay even though it is focused a substantial distance 
from its central axis. In fact, this anay is capable of focusing (o 
about 35^ from the central axis of the array without producing 
sigfiincant grating lobes. 
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So far. all of (he acoustical field measurements presented 
hav/e been in a plane parallel to the length of the array and 
normal to the surTace of the array. Ideally, (he acoustic&l field 
would taper oft toward the edges of the array but be more or 
less uniform around the arc of the array. Unfonunatcly. this 
ideal docs not match reality, as can be seen in Fig. 7, The 
acoustic field was measured as a function of rotauon angle 
along a line by fixing the position of the thermocouple m 
the focal region and roiaiipg the array about its axis. The 
acoustical intensity is not at all unifotm and. in fact, vanes 
as much as 50% before tapering off at the edges. All of the 
2-D acoustical field plots shown here were made at the 0*^ 
rotation angle where the intensity is only about 50% of the 
peak. Measurements made, but not shown here, show that this 
fluciuaiion m pressure amplitude due to the rotation angle only 
ahers the peak mtcnsity and docs not effect the overall shape 
of the acoustic field. 

Fig. Staj and (b) shows the acoustical field plots produced 
using the two multiple focus techniques: pseudo-inverse and 
split focusing. Below them, in Fig. S(c) and (d). arc thcii 
respective simulation results. Both techniques were used lo 
synthesize foci 40 mm from the array surface and 10 mm 
on either side of the central axis of the array (20 mm apan) 
using the 64-elemeni array. Notice that the pseudo-inverse 
technique produces shaiper foci than those produced by the 
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Split focusing method. This occurs because ihc pseudo-inverse 
lechnique utilizes the encirc length of che array allowing ii to 
produce sharper foci than ihe split focusing technique which, 
since it divides the array into two subarrays, effectively uses 
array lengths only half as long as (he actual length of the array. 

CM Virro Kidney Exptrimtnts 

Fig* 9(a) shows the icmperaiure r»?e vcr^u< time alons a 
fixed seven-sensor thermocouple probe located perpendicular 
to the array, The distances marked denote the distance of 
the thermocouple from the edge of the kidney neare.« to 
the array. The array, with 32 active elements, was focused 
along its central axis and 10 mm into the kidney. (Thirty- 
two active elements were u.«d since the array length in this 
conneuraiion could adequately cover the kidney.) The kidney 
was penuscd at a rate of 2.9 kg m'"* 5- ^ Note that the kidney 
is only 35 mm thick and thus, the temperatures close to the 
surfaces of the kidney are dominated by cooling caused by 



the circulating water in the water bath. All of the following 
results were obtained in the center of the kidney where the 
focus was located and the surface effects were smallest. The 
temperature rises achieved in the middle of the kidney indicate 
that therapeutic temperatures can be achieved with a staitonar>' 
single focus at realistic perfusion levels. Typical cross-section 
profile^ measured in the middle of the kidney with the pull- 
back lechnique tnr a center focus and 15 mm on either xide 
are shown in Fig. 9(b). 

The three different techniques for heating larger volumes: 
pseudo-inverse, split focusing, and scanned single foci arc 
compared in Fig. I Of a) and (b). The two multiple focusing 
techniques were used to produce foci 40 mm from the surface 
of the array and 10 mm on cither side of the array s central 
a^is, while the scanning routine produced a scries of single 
foci 40 mm from the array surface with a total span of 
20 mm. The thermocouples were located 8 mm and 13 
mm from the surface of the kidney nearest to chc array 
OS mm and Al mm from the surface of the array). The 
scanning technique produced a narrower profile than the two 
multiple focusing techniques and the distribution lacks the 
temperature drop between foci. The two multiple focusing 
techniques produce vinually identical profiles, although the 
pseudo- inverse technique does not produce as large of a 
temperature rise on (he deeper thermocouple (beyond the 
focus I as does the splii focusing technique. This t$ caused b> 
the mcreased sharpness of the foci produced by the pseudo- 
inverse lechmque. as previously explained. 



V. DtSCL'SSIOK AND SUMMARY 

.\n tntracavitar)- ultrasound phased array composed of half- 
cylinder transducer clemems has been constructed for inducing 
hypenhermia in the prostate via the rectum. A 64-channel 
amplifier system has also been designed and constructed to 
drive the phased array. As was Shown by the acoustical field 
plots, the array is capable of producing focused fields as well as 
fields containing more than one focus. The array is currentK 
capable of producing a t2°C temperature rise in a perfused 
phantom using a stationary focus, and smaller temperature 
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increases were achieved in larger volumes bv scanning the 
focus or by creating multiple suiionary foci. This was achieved 
with flows 10 the kidney ihai should simulate the relatively low 
perfusion in prosiaic (201, (21). 

What has ycc to be conclusively shown is such an orrav's 
ability to heat tissues to therapeutic tempcraiurcs in a regu- 
lar and pnsdiciable manner. The peak temperatures achieved 
withiii the perfused phantoms tended to vary considerably from 
experimcm to experiment, making direci compansons between 
ihc absolute icmperacure profiles difficuir. Therefore, the tem- 
perature profiles thai compare various focusins lechniqucs 
were normalized to better illustrate the overajl differences 
u) the shape of the icmperaiurc profiles. The fluctuations in 
temperature were caused by a variety of problems, including 
morphological differences in kidneys and rhe location of the 
thermocouples within the kidneys, the inefficiency of the array, 
and the open loop manner iti which the array was oocrated. 

The clecmcal efficiency of the anay would usuaJly drop 
considerably during the first experiment, and durin? each sub- 
sequent experimcm due to chaiiges in the eleciricai imDed^nce 
of the array elements. T^o pnmaiy factors were responsible 
for the observed changes, in the electrical impedance of the 
array elements: The first was a thermally induced impedance 
dnn caused by the array self-heatmg durine sonication. The 
second impedance shift .thai occurred from experiment to 
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expcnmcnt because the impedance of the array elements did 
nor return to ihttr anginal impedance after cooline. The 
impedance shifts made it difficult to accurately conu-ol the 
radiated power delivered to the kidney. Other problems with 
the array included arcing between elements and water seepage 
behind the anay ihat reduced the array efficiency in subsequent 
cxpenmenis until the anray was repaired and rematchcd. NVith 
a more caicfully consmjcied array, most of the power limiting 
problems experienced can be avoided. 

The acoustical field plots were in good agreement with the 
theory, and the techniques for healing larger volumes wert 
all functional, though the multiple focusing techniques were 
the most effccnve. Overall, phased arrays show considerable 
potcniial for improvement over currently used intracavitary 
ultrasound hypcnhermia system. 

While phased airavs allow signincaniJy more control over 
the acoustical field, ihe current design using half-cylinder 
radiators siill lacks conu-ol in the angular direction (around the 
arc of the array). Additionally, since the cylindrical radiators 
do not have uniform angular intensities, the angular heating 
pattern is somewhat degraded, though thennal conducuon will 
probably smooth the resulting lempcracure distribution. Similar 
fluctuations in the angular field distributions have been shown 
for other cylindrical transducers (2). (111. Radial control 
could be achieved by dividing the half-cylinder elements into 
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Ti%. 9, Tcmpertiure proftle mcitured in a perfuied kidney using 32 ictive 
elemenu focused 30 mm from ihe array surface «*-hile (he kidnev uai pcrfuied 
It me file of 2.9 kg s*"' l3 mL mm" ' flow). The amy was poiHioned 
)& mm tiom (he kidnev surr'acc. The themocoupte Idcatigni indicate the 
distance from the siufacc of the kioney. (a) The aital disifibuiion measured 
12 tnm fmm the surface of the kidiicy using a pull-back iheimocouple <to) in 
a leparaK cxpenmeni. 

pie-shaped subclcmenis and driving each subelemcnt indepen- 
dently. However, ihis would not only dramatically increase 
(he number of amplifiers, but the resulting subelements would 
have viaually no inner electrode due to the wall thickness 
of the clement. This would seriously degrade ihe already low 
acoustical efAciency of the cylindrical radiators. 

A possibly better array design would utilize planer array 
elements mounted on a rotating platform. By tilling the array 
back and forth, the same volume could be treated as with 
the cylindrical radiaiors. By controlling the acoustic neld and 
output power ai a function of tilt angle. 2-D control over 
the power deposition field csn be achieved wuhoui requiring 
additional ampliners. An additional improvement on the array 
design would be to increase the operating frequency. While the 
500-kHz operating frequency allowed larger elemeni sizes and 
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Fi| 10 Norrnaiiicd icmperaiurE/po^'er profiles measureo ir. i perfused 
ktdrtey iff muiiipic locu^inp umh^ tnc pseudo-inverse and tcix focusing 
icchniquc^. at m-cii 35 5tnek*focus ^canning; The multiple loo ^ tte created 
mm occp. IC» fT\m on Cither side of the ceniral axis whiit mc scanning 
vi\ done >^iin i ma^Ic focu^ 40 mm deep, iind scanned 10 mm on either 
vide o( the ernirat a^ia Tne fneaKun:mcnis were made t mm la) and IS mm 
(bl ffom (he surface oi' the kidney with the amy postiioned 21 mm from the 
kidney su.'face 

minimized gratine lobe formation, higher frequencici would 
increase the power absorption in tissues. Since most tumor 
sites for which this applicator would be used occur near 
iha caviiy uall. the deep penetration allowed by the 500- 
kHz operating frequency is not necessan . T^e difficulty in 
increasing (he operating frequency is that the elemeni size and 
ccnier-io-cenier eicmem spacing would need to decrease as 
the frequency mcreased. making Array design more difficult. 

As a conciu.<iion. the iniracavitar>'. elecirically focused array 
did demonstrate, ai practical depths for prostate ueatments 
t about 30 mm), the feasibility of using phased arrays for 
imracavuary hyperthermia purposes. With a more carefully 
constructed array, most of the power limiting problems ex- 
perienced can be avoided. The acoustical field plots were 
in good agreemeni wiih the theory, and the techniques for 
heating larger volume.^ were all functional, though the multiple 
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focusing techniques were the most cffcciivc. Overall, phased 
arrays show considerable improvement over ihe currently used 
intracavuaiy ultrasound hypenhermia system (I). However, it 
does not offer power control in the angular direction (around 
the axis of the applicator) and thus, new designs allowing 3-D 
conirol of the power field should be developed before clinical 
testing. 
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Cdntrbl of the necrosed tissue volume during noninvasive uiirasound' 
surgery using a IB-element phased array 
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Focused high-power ultrasound beams are well suited for noninvasive local desiruciion of deep 
target volumes. In order to avoid cavitation and to utilize only thermal tissue damage, .high fre- 
quencies { 1 -5 MHz) are used in ultrasonic surgery. However, the focal spots generated by sharply 
focused transducers become so small that only small tumors can be created in a reasonable time. 
Phased array ultrasound transducers can be employed to electronically scan a focal spot or lo 
produce multiple foci in the desired region lo increase the treated volume. In this article, ihcoreiical 
and e^pcnmental studies of spherically curved square-element phased arrays for use in ultrasonic 
surgery were performed. The simulation results were compared with cxpenmcntai results from a 
16-element array. It was shown that the phased array could conuol the necrosed tissue volume by 
using closely spaced multiple foci. The pha5ed array can also be used to enlarge a necrosed tissue 
volume in only one direction at a time. i.e. lateral or longitudinal. The spherically curved 16 
square -clcmcni phased array can proOucc useful results by varying ihc phase and amplitude scning. 
Four focal points can be easily generated with a distance of two or four wavelengths btisveen the 
two closest peaks. The maximum necrosed tissue volume generated by the array can be up to 
sixteen times the volume induced by a similar spherical transducer. Therefore the treatment time 
could be reduced compared with single transducer treatment. 



Key words: phased array, ultrasonic surgery 
I. IN7R0DUCTI0N 

In ultrasound surger>\ to avoid cavitation and have a sharp 
botmdary between the tumor and normal tissue, hjgh fre- 
quency focused transducers have to be used.' The focal spots 
generated by sharply focused ultrasound transducers are 
small when compared with the diameter of many tumors. 
Obviously this type of focus is not efficient to treat large 
tumors, which would require a large number of exposures to 
eover the whole target volume. However, in order to accu- 
ra»«iy treat the target volume close to critical structures, 
s . focal spots may be required during pan of the therapy. 
Thus, controllable focal spot size is required. The most at- 
tractive method to obtain control over the focal .^pot size is to 
use phased arrays. 

Phased array applicators were introduced to ultrasound 
hypcnhermia cancer therapy in the early 19fi0's. During the 
past decade, many efforts have been made to mvesiisate the 
advantages of phased arrays in hypenhcnmia, and several 
phased array applicators have been developed. Phased array 
applicators can be divided into the -following categories: an- 
nular or concentnc*ring arrays;^'* stacked linear-phased 
arrays/ secior-vonex arrays,' tapered linear-phased jnays.^ 
cylindrical -section arrays.^ and square-element spherical- 
section arrays. ' It was shown that ultrasound phased arrays 
can provide good control over* the heating pattern with flex- 
ibility .in moving the focus and producing multiple foci. Pre- 
vious research in ultrasound phased arrays has been concen- 
trated in hyperthermia cancer therapy. Relatively low 
frequencies (about 0.5 MHz) and small element size tcom- 
parcd to the wavelength) were employed in these studies. 
The studies on square-element sphciical-scction arrays em- 
phasized using a large number of elements, up to a few hun- 



dred. The main disadvantage associated with usins a larae 
number of elements is that the same number of amplifiers 
and electronic circuits are also required. 

The motivation of this anicle was to demonstrate that a 
phased array system u'jth a small number of elements can be 
utilized 10 provide control over the focal spot size for surgi- 
cal purposes. The excitation signals of the elements vary in 
both amplitude and phase. Hence the major usk in uulizing 
phased arrays is to determine the phases and ampiitudes 
needed to produce desired ultrasound fields. The inverse 
technique introduced by Ebbini and Cain' can be used for 
direct synthesis of ultrasound fields with multiple foci. Sev- 
eral amplitude and phase sellings ^ere calculated for differ- 
ent sets of selected control points. The limitation of the num* 
ber of utilizable conu^ol points was also investigated for a 
given array. Several simplified amplitude and phase sentngs 
based on the calculated amplitudes and phases were em- 
ployed for ultrasound field calculations. These driving signal 
sets can be utiiired. when different focal spot sizes arc re- 
quired for the array proposed here. The transient bioheat 
transfer equation was employed to estimate the temperature 
elevation due to the ultrasound power deposition. Then the 
necrosed tissue volume was predicted by the isothermal dose 
volume. Computer programs were also used to do a paramet- 
ric study to invesugate the influence of the dimensions and 
frequency of the array on the neCrO&ed ftssu£ volume. 

II. MATERIALS AND METHODS 

A. Square-element array with spherical surface 

The phased array was made from a spherically curved 
transducer which was divided into small square elements. 
The whole array was airbacked. It was constructed from a 
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FiC, I. Confi5ur4noo of ;r.f »quare-eiemcm jphcmi'.iv curved pn.vicO imy 
with lu eoordinaif tyttcm. 



Where y = p is the densiiy.of medium, e is the tpeed 
of sound, * is {he real wave number, a, = is (he 

complex excitation sousce of the nth element with amplitude 

and phase 6^ . is the disunce froir a point ix„ 
on the fiih element to the field pomt of interest ix„ »y^,£„). 
5^ is the area of the nth element. 

In uUrasonic surgery, waicr is used as z coupling medium 
between the transducer and tissue interface. Waier is consid- 
ered a nonatienuatmg/medium. and tissue a weakly anenuat- 
inff medium. Tne effects of anenuation can be saiisfactoniy 
described by replacing e"*'^'*** with ^ m expression 
(M. where k^^k-ja is the complex wave number wiih at- 
tenuation eoefneicni a.'* 

The power depostuon for the desired volume ts given by 

|pmU.y.ci|- 

GlJr.\.^J = M^^. z . (2j 



spherically curved PZT4 bowl, having a diameter of 100 
mm. a radius of curvature of 130 mm, and a frequency of 1 .4 
MHz. The bowl was cut into 16 elements, each with a length 
of 20 mm per side, a 0.3-mm space between the elements 
was filled wiih silicone rubber for electrical and mechanical 
isolation. Each of the elements was connected to an LC 
matchinff circuit to match the impedance ro 50 O and 0*. Tlie 
array was driven \Mth a custom made 16 channel ampliner 
(Labthermics. Champaign. Illinois). The phase and amph- 
cude were controlled by RF signals feeding the amplifiers. 
l^t driving signals ^ere generated by an ni-house manufac 
tured digital circuit.* The amplitiide and phase of each mpui 
signal were digitaii\ controlled with 3-bit resolution in am- 
plitude and 4-bii resolution m phase. The number of ck- 
ments in each row ano each column were the same so ihaf 
the whole applicator uas a square shaped focused transducer. 
Tne configuration of the phased array for the cxpenments m 
shown in Fib. 1 . 



B. Ultrasound field measurements 

The relative pressure amplitude squared distributions 
were measured in as gassed water using a needle hydrophone 
(active spot size 1 mmi scanned across the focal region. The 
needle hydrophone was moved by stepper motors, typically 
with 0.1 -mm steps siross the beam. The total acou.stic power 
was measured usin; a radiation force technique. 



C. Computation ot ultrasonic fields • 

Consider an uiirafcund phased array wiih N ciemenii. 
employed to produce ;tn ultrasonic field in a nonattenuaiinf 
medium. Assume tne coordinate -.system is defined as in Fig. 
1. According to (he theory of ultrasonic radiators developed 
by O Neil'" and the principle of superposition of acoustic 
pres.^ure. the prepare field due to this transducer can be 
evaluated by ihc'Ra\leich-Sommerfcld integral The inie- 
cral IS 



where /t,^, is the absorption coefficient and Z = pc is the 
impedance of the medium. 



0. Inverse technique 

The mvcric lechniquc can be used lo calculate the ampli- 
tude anc pna^^e sellings from selectee control points wnere 
the desireo pressure values are given. Defining 

JOCK I c ' — • 



tncn for .^1 field points, expression ili can c: written :r 
matrix notation as 



f4 



i P\ 

I • ' 

I Pw 



U t 



H 



U K 



The elements of matrix H are evaluated by numerical intC' 
gration using expression (3). Equation (4'i has \vto imponani 
features. Fir&t. if the H mamx is caleuiaied ar.c saved for a 
desired field, then the pressure fie!d can be evaluated using 
Eq. instead ofEq. (1 ) for a given excitation source manri.x 
U. Compuiaiion time can be saved by using iriis method 
when vQr)ing the cxciioiion source over the same caiculauor 
volun:. Second, for a ocsired pressure field pziicm, i.e.. for 
3 given matnx P. :hc excitation source U car. c: calculated 

b> 

wnere i^ the p.seudoinvme mairix of H. L*5ually the loul 
number of control points arc much le.«s than the loiai number 
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Table I. $«lc^«d comrol puinm uwO fur invert calcubwonss for 16 »quarf 
*' cleincnc spficricsl curved pha>cd array. The eonirtil polnm ire loCitieJ «i 
* :«*129 nun plane. Note ih» uniu .» and T -Jit millimeter. 
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of exciuijon sources in ihis case. The W mairix can be 
foi»*<l by uiilizins single value decorriposinon. 

.t control points were selected in a plane ai a desired 
focal position. All the points were evenly distributed on a 
circle. Since the pressure at a control point is a complex 
number, both the magnitude and the phase are required to 
perform the inverse calculations. In selecting control points 
in this manner, it is natural to make the magnitudes the same 
for all the points, and the phases either in phase or out of 
phase so that the phase rotates around the central axis. The 
latter setup gives destructive interference on the central axis 
thus, eliminating potcmial hot spots on the axis.' The sc- 
Icacd control points used in ihc inverse calculations ar: 
given in Table 1. 

E. Thermal modsling 

An approximate temperature re»pon.^e to the power depo- 
sition vvas predicted by the transieni bioheat transfer equa- 
tion. The differential equation is 



where T is the temperature Vimei. . XlfbSxailoA ('ix^/.z!!.. 
'pf is the density of the tissue, c, is the specific heal of the 
tissue, li, is the thermal conductivity of the tissue, u* is the 
blood perfusion rai« . is th«. specific heal of .the bloci. T„ 
is the arterial blood temperature, and Q{x.\\z) is the acous- 
tic power deposition rate per unit macs. The ihermal respon.^e 
was simulated in a homogeneous medium. A surface tem- 
peraiure of 37 'C on (he cube and an initial temperature of 
37 *C inside the cube were used as boundary and initial con- 
ditions for all of the compuiaiions. A numerical finite differ- 
ence method was employed to soIvIe Eq. (5). Previous studies 
have shown :hat there arc sevcril parameters which affect 
the temperature elevation.'"^ Temperature elevation has been 
shown (0 be almost mdependcm of the blood perfusion rate 
for ihon ultrasound pulses.'^ for this reason, ultrasound 
pulse durations of 1. 5. ard 10 s were used in this study. The 
maximum temperature reached in all simulations was kept 
the same (80 *C) by adiustine the input power. 

F. Thermal doso calculation end necrosed tissue 
volume estimation 

The thermal dose calculation was based on the technique 
suggested by Saparcto and Dewey. '-^ U$ing this technique, 
the accumulated th^ermal dose was calculated at a reference 
lempcraturc by numerical imearaiion under different tern* 
pcraturc profiles. The ihcrmal dose. i.e.. equivalent time, ai 
the reference temperature can be evaluated by 

Dosei/r )= dt-' .2 /^^-r^A^r. 



where T",,/ is the reference temperature, /f., 

is the final nme. ^: js a smyll time interval. 7"^, is the avct- 

asc tcmocrature dunns time ^i. and is a paiamcicr giver 

by 

'0.3. if rn)&43 
\0.2^. otherwise 

The necrosed tissue volume is estimated by the. isothermal 
dose volume surrounded for 240 mm at the reference tem- 
perature of 43 °C. This technique has been found to be i 
reasonable model for predicting tissue necroses induced by a 
spherically curved transducer.''^ 

A two layered medium water-tissue was assumed in the 
simulations. The speed of sound and the density were 1500 
m>s and 998 kg/m' , respectively, for both media. The attenu- 
ation coefficient of the tissue was assumed \0 be 10 Np/m' 
MHz. The thermal properties of the iis.sUc are given in Tabic 
II. 

ill. RESULTS 

A. Comparison of expenmental results with 
simulations 

It was nccessarv- to venfy the numencal model using the 
16 square-element phased array with a spherical surface be- 
fore relying on the simulations. The cxperimcniai and simu- 
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lated profiles of ihe pressure amplitude squared at the acous- 
tic focus are displayed in Fig. 2. All (he data were 
nonnalized to 1 by dividing by the maximum pressure am- 
piimde squared for each cun-c. Figure 2(a) shows ihc resuiis 
for the case a*iih uniform excitaiion sources. Tnere was sood 
agreement bec\i-een the simuiatiorvs and expcnmcmai resuiis 
for ihc mam beam. When the array was excited by various 
phase setanes. the simulations reasonably predicted the side- 
lobe distributions found in the experiments [Figs. 2(b)-2(d)) 
The maximum power output measured was at least 14 W for 
uch transducer element. 

B. Simulation re9utts 

4 

The phased array ultrasound transducer used in the ex- 
periments was simulated first. A calculation volume of 30 
X3QX90 mm^ wai used for this phased array. The axial dis- 
unce from the center of the transducer to the water-tissue 
interface was 70 mm. Tne inverse technique results arc 
shown nrst. The amplitude and phase settings calculated b> 
the mvcrsc technique corresponding to selected control 
points arc given in Tabic III (only the first four are shown. 
Tne contour plots of power deposition calculated for these 
anay settings are dispiaycd in Fig. 3. For four control points 
only one focus was produced., which was centered on the 
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axis when (he pressure at all of the conirol points were in 
phase [Figs. Mil and 3(b)]. Four focal spots were obtained if 
the pressure at the coniroi points had 90.0" phase shifts from 
one point to another [Figs. 3(cl and 3(d)] Comparing Figs. 
3(c) and 3(dt with Figs. 3iei and 3(0, tha peax to peak dis- 
tance was increased when the control points on a 1.25*mm 
radius circle increased to a 3.25-mm radius circle. However, 
if tne cor.tioi points were on a S.25'mm radius circle, the 
tocal spots were at the same location as Figs. 3fci and 3(d). 
Tne ncid oisinbunons mot shown herci were very similar to 
Figs J CI and 3<d) For Six control points. s:x focal spois 
were obtainec at control locations with different panenis 
[Figs. 3igj and 3(h)]. sno two additional focal spots were zt 
the center. When the number of control points increased to 
eizhi. the power deposition pattern had the maximum pres- 
sure ampliiiide ;n irohi of the focal plane w:ui strong peaks 
at the watcr-iissue inienacc [Figs. 3(i) and 3(j;]. Tncsc re- 
sults were used to select simplified amplitude and phase set- 
tings wh:ch ^>^cre then used in necrosed tissue voiume simu- 
lations (Table IVj. 

The isothermal doses for tissue necrosis tor live amplitude 
and phase seuinn arc aispiayed in Fig. 4. Let us consider the 
ultrasound puhe duration lO-s case (solid line'i. For the uni- 
form excitation case, the focus was located at a depth of 59 
mm rrom ih: interface. Based on the 240-min isoiherma: 
dose [Figs. 4iai and 4rb)]. the shape of the necrosed tissue 
volume was ciose to an ellipsoid. Tne calculated volume was 
about 85 mm*, measunng 3 mm laterally and 18 mm longi- 
tudinally. Isothermal doses ari: also shown in the same hg- 
urcs for the ultrasound pulse durations of 5 and 1 $. where 
the input power was adjusted so that the maximum tempera- 
ture was Icepi at SO 'C. For case 1 of Table IV. ihe necrosed 
nssuc length (Fig -J'ci] was 42 mm. and the width [Fig. 
4(di) was 3.4 mm. The length was more than double and the 
widih was slightly enlarged compared to the uniform excita- 
tion, case. Tht computed ncc''0S''<* tissue volume was ibou: 
339 mm'\ Ficures 4ic) and 4(0 snow case II ot Table IV. for 
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FW. 3. Coniour plou of power depojinon tor cilculaied wnphiudc and pnase itiiingj. Th' -)ntrr> poimi usca to cilcuUie unpliiudes and phues from Table 
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Table IV. Tht simplified phue ind amplitudt scttmis foi the I6 squuv- 
. element sphefieally curved phiSM vny Ufctf in the necrosed Uiiuc volone 
caJculitioAS. 
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which ihe necrosed tissue length and width were 20.8 and 
5. 1 mm. rcspccuvdy. Tne length was slightly longer and the 
width was almost double compared wjih the unifonn cxcita- 
lion case. Notice that this case is equivalent to the siiuation 
in which the phased array had only four elements with sides 
of 39.4 ram. The calculated necrosed tissue volume was 
about 377 mm'. Figures 4(g) and 4(h) mdicatc that there 
were four focal points in the power deposition for case III oi 
Table IV. The shortest peak to peak distance was abour 4\. 
where X is the waveiengih. Tne whole voiumc u a$ ihe suit.- 
mation of the four tndividuaj smali necrosed tissue volumcf 
which were produced by the four focai pomis. Each smai; 
volume was about the same .-ize as [he voiumc m the uni- 
form excitation case for the I- and^i-s sonicaiions. For the 
10-s jonicaiion a unufed volume was created wnh the total 
volume of 810 mm*. For the amplitude and phase seiunj of 
cast rv in Tabl: IV. the power deposition tnoi shown here; 
indicated that mere were four strong focai points tdepth 54 
mm I surrounded by tour small focal points. The shortest dis- 
tance between the strong peair« about 2X. SimuhaneouA 
focusing at these points was penonncd to enlarge ine heated 
volume compared to the uniform excitation ca.^c. The calcu- 
lated isothermal dose (Figs. 4ii) and A(\)] from the tempera- 
ture distribution had a maximum length of ?8.5 mm and i 
maximum width of 7.2 mm. and the necrosed tissue volume 
was close to 1 .4 cm". 

To illustrate the effect of frequency on necrosed tissut 
volume produciion, chc isothermal do5es for case IV of Table 
IV wiih an ultrasound pulse duration 10 s zje sriown in Fiss. 
fia) and 5ib) Ths calculated volume vvas 60x60x90 tT)r>- 
for 0.5 NfHz. 40X40X90 mm- tor 1.0 MHi. and 30x30x90 
mm* for 1.4 and 2.0 .VlHz. The distance from i/ie transducer 
to tne interface v^as kept the same (70 mmi. while the f:t- 
Quency was vancd. The ratio of the necrosed tissue length to 
widtn A-as almost independent of frequency. T.ie iize in- 
creased as the frequency decreased. Tne shapes were similar 
for frequencies of I ic 2 MHz. 

40 
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To understand the effect of radius of curvature on the 
necrosed tissue volume, the isothermal doses for case IV of 
Table |V. with an ultrasound pulse 'duration of 10 s and op- 
erating frequency of 1.4 MHz aic displayed in Figs. 5(c) and 
5(d). The distance from the transducer to the interface was 
varied according to the radius of curvature to maintain a 
constanc focal depth. Tne necrosed tissue volume increased 
rapidly along the longitudinal direction toward the interface 
as the radius of curvature increased. However, the diameter 
of the necrosed tissue volume increased only slightly in the 
focal plane. 

To 4how sne capability of the phased array to move the 
focus, contour plots of the power deposition for tne maxi- 
mum displacement arc displayed in Fig. 6. Tne simulation 
model ^as tne same as the one used to generate Fig. 3, The 
phase distributions were obtained based ~on direct cajcula- 

lions. J e.. phQst c». = 2Tr</„/X. /i 1 .2 .V. where is ihc 

disunce from the center of the element to a desired focus. 
Figure 6ta) shows the results for (he case with uniform am- 
piituoe and pnase setting, the focus could be moved 7 mm 
closer to me interface [Fig. 6(b)). or 7 mm deeper [Fig. 6(c)] 
compared wni) Fig. 6/a). Figure 6(d) shows the phased array 
focus shit led I 3 mm off the central axis. When larger dis- 
placemer.ii were attempted, the phase increment between- ad- 
jacent eiemenTs exceeded 7r/2. Distributions generated by 
those onascs are no longer single, strongly focused ulua- 
sound ne.z! rr.crcforc these were the ma.ximum displace 
menis acnievabie with this array in the sense of keeping the 
fields jn ihc rs sion of fccus similar lo ihcie fcr the uniform 
excitation case 

Finally. r\;, 7 showj ihe input power requirements for the 
uniiom c\c!;at:on and four amplitude and phase 5enin?s 
piven ifi Tas.c iV For Ultrasound pulse durations cf 5 or 10 
s. inc powder ranged from 30 to 250 W For a i-s pulse, the 
input pouer could be as high as 900 V»'. depending on ihs 
ampntudr ur.c phai: sellings. 

IV. DISCUSSION 

The exper.menial and simulation results shoMbed (hat 2 
16'ciemen: pnased array can offer significant ccmroi o%cr 
the size and shape of the necrosed ussue volume dunne u!- 
trasouna furgery. Tne experiments showed that sucn an arra) 
can be ccnsiructed anii 11 can deliver enough power for sur- 
gical purposes. The simulated and experimental ultrasound 
field distnDu:ions aereid reasonably well. The differences 
may dc due :c uneven element size and power output m the , 
phased arav vvcii as mca5urerT>er\i errors. Since the peak 
to peak dis:ancc was only a few millimeters, the ultrasound 
dctccic: T.jv contnbutc some measurement errors. The simu- 
lation mode; accurately predicted the locations of sidelobes 
and main beams except that 11 over predicted :he magnitude 
of the sidclob£5 When adiacent elements of a ohased array 
are ver\- ciose to each otner. the mutual eoupiing between the 
elements would have an efi'ect on the ultrasound neld. The 
difference oetween tne sidelobes in the expenments and the 
simulations may be due to the lack of such crosstalk in the 
calculations 

The ouove resuiu proved (hat the large element spn^n* 
cally curved phased array can enlarge the necrosed (issue 
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volume. The phased array can also enlarge the necrosed us- 
suft volume in only one dircciion ai a time, if desired. Ji is 
important to be able to control the focal spot size so that 2 
large tumor could be treated in a reasonable time. The con> 
itniciion of the whole system was relatively simple due to 
the small number of elements. The 16-elcmcm phased array 
can aiso shift the focus off the central axis or move the focus 
along the central axis. However the distance is limited to 
^1.5 mm laieraUy and 14 mm in the axial direction. In mov- 
ing the focus along the central axis, the array is similar to 2 
concentric-rins array with two rings. Theoretically, ihc maxi- 
mum phase increment between adiacent elements is r: 
Therefore, the maximum possible phase difference between 
the smallest and largest phaiscs is it when moving the focus 
along the central axis. From geometric considerations, this 
phase difference produces displacements along the central 
axis of up to 23 mm (10 mm closer, 13 mm deeper). In 
shifting the focus sideways, it is simitar to a cyhndrical- 
secuon array with four elements. The maximum possible 
pha.se difference between the smallest and largest phases is 
3 IT for shining.ihe focus sideways. Geometrically, ihis phase 
difference Shifts the focus r3.4 mm laieraily. However, the 
simulations snowed (hat the phase increment between the 
adjacent elements shou' ' be less th>n it/2 (o generate n 
single strongly focused ultrasound field. 



The 16-eiemem phased array can generate four focal 
points with a pcaK to p^ax distance as short as (Wo wave- 
lengths The maximum distance bef«.'ccn the closest peaics is 
limited 10 acoui four waveleneihs. When the four control 
points were on a circle of radius 5.25 mm, the phase incrc- 
mcnt beiween jdiacent elements exceeded tt for this array 
Therefore, it !s physically impossible to use the phase and 
amplitude setting obtained by inverse techniques to produce 
the four tocai spots ai :'ne comrol points. Although Eq (4) 
represents an unoerdetermind system, it produced a solution 
at the control noints When multiple foci are separated by 
distances larger than 4.6 mm. the necrosed tissue volume 
becomes a feu/ individual smaller volumes, which arc not 
united for short ultrasound oulse duranons 

Bv decreasing the frequency, the necrosed tissue volume 
can be eniarccd- because the focal Spot increaies due to the 
increased wavelength. But the ratio of the necrosed tissue 
length to the width was kept almost the same. As the radios 
of curvature is increased, the necrosed tissue volume can also 
be enlarged. The necrosed tissue voiume is increased mainly 
in the axiai direcnon. This agrees with previous experience 
using single element spnerically curbed iransoucers.'' 

Phased atrays require more input power than similar 
'tnele focuseo transdiicrrs to reacn the same temperature 
level for the same ultrasonic pulse duration due to increased 
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focal spol size. The test array measurements showed that 
Uiercquired power levels can be generated in practice ai least 
^- and 10-s exposures and thus, the proposed technique is 
It^sible. 

Irt order to obiAtn the desired ijeld patieni. several field 
control points are necessary 10 perform the inverse calcula- 
tions. However, the number of uiiliiablc field control points 
is limited by the number of phased array elements. Hence the 
number of field panems that produce significantly different 
shapes or sizes (]f the measured volume is limiied. .Although 
only few amplitude and phase settings were presented in this 
study, (he results provided information which could aid in 
planning ultrasonic s.ursery using a phased array system with 
a small number of elenents. These resulis have not been 
optimized: nevertheless ihey illustrate the feasibility and 
range of focal size and shape obtainable with a 16* element 
array. The calculations and experiments can be used 10 obtain 
a set of different focal spots for a given array. This limited 
number of focal spots can then be utilized 10 optimally cover 
the target volume in minimum time. The mi^imi^alion of the 
duration of the treatment is an imponant factor for control- 
ling the cost of the procedure, especially when MR I is used 
to guide and monitor the treatment.'^ If more control over 
the field panem is required then an anay with a larger num- 
ber of elements is required. Such an array could also be used 



I 

c 



900 










800 


: 10 s 






z. 


700 


g 1 s 






t, • 
*, t 
• • 


600 








• « 

* ' * 


500 






* 


• « 

♦ • 








1 • 

* 




400 




* • 

* • 

•« • 

m • 


• 


< . . 


300 


• 




• • 

* 


» » 

m 

• « 


200 
100 


• 

• ■ 

* ft 




* 

• _* ■ 

• •• * 
• 


t «« » ' i 






■ 'Cv 




0 












CBse ^ case 2 


cases 


case 4 


csisa S 



Ftc. The input po«»'er requucmsm to reach ihc &ame tcirperature le^ti 11 
the end of uliruound puhci for ^anous caaec^Ca^c I u for the unifam 
eaeiution caic. cast 2-j arc corrcspondtn; 10 the implituda and pha.«s 
seningt of cast l-|V in Table IV. 



to compensate for the effect of tissue nonuniformitics which 
may distort the focal pauern.'' *' 
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An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 

a plurality of transducers, 

an excitation source for driving the plurality of transducers to generate and transmit 
ultrasound through the skull to induce cavitation at least at a selected region of the brain, 

the excitation source driving at least selected transducers at differing phases with respect 
to one another. 

i 

An. apparatus according to claim 1, wherein 

the excitation source drives each of the selected transducers at a phase that compensates 
for a phase shift effected by the skull in the ultrasound generated and transmitted to the 
selected region by that transducer, 

so that the ultrasound generated by the selected transducers arrive substantially in phase 
with one another at the selected region. 

An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
transducers arrive at phases within 90^ of one another at the selected region. 
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An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
transducers arrive at phases within 45** of one another at the selected region. 

An apparatus according to claim 2, wherein the excitation source drives each of the 
selected transducers at a phase so that the ultrasound generated by the selected 
transducers arrive at phases within 20° of one another at the selected region. 

An apparatus according to claim 1, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.01 
MHz to 10 MHz. 

An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0. 1 
MHzto2MHz. . 

An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

An apparatus according to claim 8, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

An apparatus according to claim 9, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 
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an ultrasound transducer, 

an excitation source for driving the transducer to generate and transmit ultrasound 
through the skull to induce cavitation at least at a selected region of the brain. 

An apparatus according to claim 1 1, wherein the excitation source drives the transducer 
to deliver ultrasound to the selected region at a frequency ranging from 20 kHz to 10 
MHz. 

An apparatus according to claim 6, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.02 
MHz to 10 MHz, 

An apparatus according to claim 13, wherein the excitation source drives the plurality of 
transducers to. deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

An apparatus according to claim 13, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

An apparatus according to claim 15, wherein the excitation source drives the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising the steps of 
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placing a plurality of transducers in the vicinity of an exterior surface of the patient's 
skull, 

driving the plurality of transducers to generate and transmit ultrasound through the skull 
to induce cavitation at least at a selected region of the brain, 

the driving step including driving at least selected transducers at differing phases with 
respect to one another. 

A method according to claim 17, wherein step (C) includes driving each of the selected 
transducers at a phase that compensates for a phase shift effected by the skull in the 
ultrasound generated and transmitted to the selected region by that transducer, such that 
the ultrasound generated by the selected transducers arrive substantially in phase with 
one another at the selected region. 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 90° of one another at the selected region. 

4 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 45*' of one another at the selected region. 

A method according to claim 18, wherein step (C) includes driving each of the selected 
transducers at' a phase so that the ultrasound generated by the selected transducers arrive 
at phases within 20° of one another at the selected region. 
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A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.02 
MHz to 10 MHz. 

A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.1 
MHzto2MHz. 

A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

A method according to. claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

A method according to claim 17, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 

placing an ultrasound transducer in a vicinity of an exterior surface of the patient's skull, 

driving the transducer to generate and transmit ultrasound through the skull to induce 
cavitation at least at a selected region of the brain. 
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A method claim 27, wherein step (B) includes driving the transducer to deliver 
ultrasound to the selected region at a frequency ranging from 0.02 MHz to 10 MHz. 

A method according to claim 28, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a frequency ranging from 0.1 
MHz to 2 MHz. 

A method according to claim 29, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region at a sonication duration ranging 
from 100 nanoseconds to 30 minutes. 

A method according to claim 30, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with continuous wave operation. 

A method according to claim 31, wherein step (B) includes driving the plurality of 
transducers to deliver ultrasound to the selected region with burst mode operation, where 
burst mode repetition varies from 0.01 Hz to 1 MHz. 

An apparatus for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising 

a plurality of transducers arranged in a two-dimensional array, and 

an excitation source for driving the plurality of transducers to generate and transmit 
ultrasound through the skull to induce cavitation at least at a selected region of the brain. 

An apparatus according to claim 33, wherein the plurality of transducers are arranged in 
an array having any of a substantially circular and a substantially annular cross-section. 
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An apparatus according to claim 34, wherein the plurality of transducers are arranged in 
an array having a substantially circular cross-section. 

An apparatus according to claim 34, wherein the plurality of transducers are at least one 
of (i) mounted in, and (ii) separated from one another by, a damping agent. 

An apparatus according to claim 36, wherein the damping agent is any of a natural or 
synthetic rubber. 

An apparatus according to claim 34, wherein the excitation source drives at least selected 
transducers at differing phases with respect to one another. 

An apparatus according to claim 1, wherein excitation source drives each of the selected 
transducers at a phase that compensates for a phase shift effected by the skull in the 
ultrasound generated and transmitted to the selected region by that transducer, so that the 
ultrasound generated by the selected transducers arrive substantially in phase with one 
another at the selected region. 

A method for delivering ultrasound through the skull to the brain of a patient, the 
apparatus comprising the steps of 

placing a plurality of transducers in the vicinity of an exterior surface of the patient's 
skull, 

I 

driving each of at least selected transducers generate and transmit ultrasound through the 
skull, and determining a phase shift effected by the skull in ultrasound generated and 
transmitted by each such transducer. 



51 



wo 98/07373 



PCT/US97/14760 



driving the plurality of transducers, together, to generate and transmit ultrasound through 
the skull to induce cavitation at least at a selected region of the brain, 

the driving step including driving at least selected the transducers at phases determined in 
• accord with step (B) so that the ultrasound generated by at least the selected transducers 
arrive substantially in phase with one another at the selected region. 

A method according to claim 40, wherein the selected region ranges from 1 nmi^ - 1 cm^ 
in volume. 
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